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Abstract 



Motivated by the concept of ideal mutual statistics, we study a multi- 
species Calogero-Sutherland model in which the interaction parameters and 
masses satisfy some specific relations. The ground state is exactly solvable if 
those relations hold, both on a circle and on a line with a simple harmonic po- 
tential. In the latter case, the one-particle densities can be obtained using a 
generalization of the Thomas-Fermi method. We calculate the second virial 
coefficients in the high temperature expansion for the pressure. We show 
that the low-energy excitations are the same as those of a Gaussian confor- 
mal field theory. Finally, we discuss similar relations between the statistics 
parameters and charges for a multispecies anyon model in a magnetic field. 



PACS Nos.: 5.30.-d, 71.45.Jp, 74.20.Kk. 



Keywords: Calogero-Sutherland model. Exclusion statistics, Anyons. 



E-mail: diptiman@cts.iisc.ernet.in 



1 



1. Introduction 



The Calogero- Sutherland (CS) model |]T], 0| consists of a system of identi- 
cal particles in one dimension which interact with each other via a two-body 
inverse square potential. The model has been studied for many years from 
several points of view. It is exactly solvable and integrable 0, |], it has 
interesting connections with random matrix theory and it provides an 
example of exclusion statistics 0. 

The concept of exclusion statistics introduced by Haldane has attracted 
considerable attention recently [7-15]. The examples of ideal exclusion statis- 
tics known so far include the Calogero- Sutherland (CS) model, spinous in the 
Haldane- Shastry model [0, and anyons residing in the lowest Landau level 
in a strong magnetic field [7-10]. The term ideal means that the only inter- 
actions between particles are statistical, and that these interactions are scale 



invariant, i.e., independent of the momentum or energy scale |]14 . 

The idea of exclusion statistics can be generalized from systems with a 
single species of indistinguishable particles to systems with several species. 
In the latter case, the concept of mutual statistics has proved to be very 
useful [10-12,14]. An example of mutual statistics is a multianyon system in 
two dimensions [§, |lOl- It seems interesting to ask whether there is a one- 
dimensional model which can realize mutual statistics . With this and other 



motivations, a multispecies CS model was examined in Ref . ||T^ ; however the 
model involved both two-body and three-body interactions. In this paper, we 
ask whether it may be possible to realize mutual statistics in a CS model with 
only two-body interactions. Somewhat unexpectedly, we discover that ideal 
mutual statistics cannot be realized unless the interaction parameters and 
masses of the different species satisfy certain relations. However we are not 
able to prove that the model does exhibit mutual statistics if those relations 
are satisfied. 

The outline of our paper is as follows. In Sec. 2, we present a general 
analysis of a model satisfying ideal mutual statistics on a line. We assume 
that the mutual statistics is described by a symmetric matrix. (This assump- 
tion is motivated by a CS model). We find that the masses and interaction 
parameters must satisfy certain relations in order that the asymptotic Bethe 
ansatz be applicable to the ground state. In Sec. 3, we consider a CS model 
which satisfies those relations. It is known that the ground state is exactly 
solvable for the model on a circle; it is also exactly solvable if the CS model 
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is placed on a line in an external simple harmonic potential, with the same 
harmonic frequency for all species 0]. We discuss a class of exactly solv- 
able excited states, including some subtleties involving the centre of mass 
excitations on a circle. In Sec. 4, we compute the one-particle densities in 
the harmonic potential problem using a generalization of the Thomas-Fermi 
method. In Sec. 5, we study the CS model at high temperature, and cal- 
culate the second virial coefficients in a high temperature expansion of the 
equation of state. In Sec. 6, we study the low-energy spectrum and eigen- 
states to leading order in 1/L. We find that the spectrum is the same as 
that of a c = 1 Gaussian conformal field theory. We derive the leading terms 
in the specific heat and pressure at low temperature. Finally, in Sec. 7, we 
consider a multispecies anyon model in an uniform magnetic field at zero 
temperature and show that it has some desirable properties if the charges 
and statistics parameters satisfy the same kinds of relations 



2. Mutual Exclusion Statistics in One Dimension 

Given a many-body quantum system with a finite number of states and 
A species of particles, let Na denote the number of particles of species a 
(a = 1, 2, y4). Let Da be the number of states which are available to the 
next a particle which may be added to the system. The -D^'s are functions 
of the numbers Na- The mutual statistics parameters Qah are defined as 

= - 9ab AN, . (1) 

b 

If the diagonal matrix element Qaa = (or 1), the a particles are bosons (or 
fermions). We will henceforth assume that all Qab > 0. In the literature, 
the off-diagonal element gab is referred to as the 'mutual' statistics between 
species a and b. 

For an ideal system with a single species such as a CS model, the param- 
eter g = gaa can be interpreted in terms of the phase shift which appears in 
a two-particle scattering [|1T|. Namely, the scattering phase shift is given by 
5{p) = 7rg6{p), where p is the relative momentum of the two particles, and 
6{p) = p/\p\ for p ^ 0. (Thus the wave function of two bosons is symmetric 
with g = 0, and that of two fermions is antisymmetric with g = 1). At zero 
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temperature, the particles occupy momenta satisfying 



P^ = Oip.-p,), (2) 

where L is the length of the system. Eq. may be derived by an asymptotic 
Bethe ansatz analysis of the CS model on a circle of length L |Ty| . Thus 
momenta from —pp to pp are occupied with a spacing of 27chg/L. If is 
the number of particles, we have 

VF L 

pp ^ 2iThg ^ 

We are eventually interested in the thermodynamic limit N,L ^ oo keeping 
the particle density p = N/ L fixed. Then 

Pf = Tchgp , (4) 

and the total energy 

E (i vr^^^^f^ ^ 3 

J-PF 27rhg 2m 6m 

We may now generalize the phase shift argument to a multispecies model. 
Let us assume that a scattering between two particles i and j produces a 
phase shift 7rgij6{vi — vj), where v = p/mis the velocity, m^mj/ (rrij + mj) is 
the reduced mass, and (fj — Vj)mimj/{mi + mj) is the relative momentum 
of the two particles. The two particles may or may not belong to the same 
species; gij = gab, = ma and mj = m^ if i and j belong to species a and 
b respectively. (It is clear that gab = gba in such a model). If we assume the 
asymptotic Bethe ansatz to be valid for the ground state of this model, we 
obtain ^ 

rriiVi = -j- 9ij Sivi - Vj) . (6) 

We assume that no two particles have the same velocity; this will be verified 
in Eq. (P) below. 

If = J2aNa, Eq. d) has A^! /UaNJ distinct solutions; this is be- 
cause the particles 1,2,...,A^ can be ordered with increasing velocities in 
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that many different ways. (Recall that different orderings of identical parti- 
cles are physicaly indistinguishable). The asymptotic Bethe ansatz ground 
state wave function will then be a superposition of that many different plane 
waves. For this idea to work, it is clearly essential that the total energy 
E = Y.if^ivfl'^ be the same for all the different solutions of (|]). We find 
that this is true if and only if 

9ij9ik 9ij9jk /_\ 

nii rrij 

for any three different particles i,j and k; they may or may not belong to 
the same species. The non-trivial solution of (|^ is given by 

Qij = a miTTij . (8) 

Eq. (|^) is one of our main results. Note that the A x A matrix Qat then has 
only one non-zero eigenvalue. 

In deriving from we have assumed that all the off-diagonal ele- 
ments gab > 0. There are other solutions of (|^) in which an off-diagonal gab 
vanishes. However in that case, (|^ implies that gac9bc = for all c. We 
can then prove, by induction in the number A, that the system consists of 
two or more disjoint groups of species, such that gab = whenever a belongs 
to one group and b belongs to another. We can therefore study each group 
separately. We will not consider such cases any further, and will assume 
henceforth that all the gab satisfy (||). 

It is rather surprising and novel that mass- dependent conditions such as 
® appear to be necessary for mutual statistics. Such conditions may be 
peculiar to models in which the two-particle scattering phase shift does not 
depend on the magnitude of Vi — Vj. We do not know of any physical system 
in which such relations might arise in a natural way. 

Using (H), we can solve (|^) explicitly. If the particles are ordered so that 
vi < V2 < ... < vn, we get 

V, = ^Y. m,0{^-J). (9) 



From Eqs. and (|), we deduce the ground state energy to be 

6L2 



Eo = [ ( E ^aNa )^ - E ^iNa ] . (lO) 
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The Fermi momentum for species a may be obtained in two different ways. 
Firstly, for an ideal system at zero temperature, the chemical potential fia = 
p\^/2ma- Further, dEo/dNa = fia- In the thermodynamic limit, the Fermi 
momentum for species a is therefore given by 

PFa = nh J2 9ab Pb , (11) 

b 

where pb = N^/L. We observe that the Fermi velocity 

Vp = ^— = Ti'ha rrib pb (12) 
ma V 

is the same for all species. The second way of deriving vp is to examine (|^); 
amongst all the solutions of that equation, the one with the largest possible 
value of Va is clearly given by (|1^). We note that the ground state pressure 
and compressibility are given by 

dL Sirha ' 

We can understand Eq. (|1|) semiclassically by identifying the number of 
states with the 'volume' of phase space. The volume of phase space occupied 
in the ground state by particles of species a is given by / {dxdp/27ih) [l+9{pa — 
ha)]/2, where ha = /2ma is the single-particle Hamiltonian. The volume 
available to the next a particle is then given by the 'remaining' amount 

Da{t) = / ^ ^[1 - d{Pa - ha)] exp(-e/i„) , (14) 

where we have introduced a cutoff e to make sense of the divergent integral; 
we will let e — > 0^ at the end of the calculation. (The final result should not 
depend on the cutoff procedure [§). From Eq. ([TT|) , the change in Fermi 
momentum is given by 

^PFa = ^ J2 9ab AiVfe . (15) 



P = - 
K = — 
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Combining Eqs. ([T^-[T5|), we find that ADq = lim^^Q ADa{e) is given by Eq. 
(i- 

In the entire discussion above, we have imphcitly assumed the system to 
be homogeneous, that is, the densities pa do not depend on the position x. 
This is a non-trivial assumption for the ground state of a multispecies model. 
For instance, consider a two-species model in which gi2 is much larger than 
both ^fii and g22- Then it is intuitively clear that the ground state will exhibit 
phase separation, i.e., the particles belonging to species 1 and 2 will prefer to 
reside in two different regions of apace. One would then use Eq. @ in the 
two regions separately, where pa is the local density of the species appropriate 
to a particular region; it would be incorrect to use the multispecies Eq. (pJ]) 
where pa is the average density defined over the entire system. 

One can derive the leading term in Eq. ([T0|) in a different way using 
the result (|^) for the case of a single species. For a CS model, the particles 
have such a singular two-particle interaction that no two of them can cross 
each other. (This will be explained in more detail in the next section). The 
Hilbert space thus breaks up into a number of different sectors corresponding 
to different orderings of the particles on the line. It then seems reasonable to 
demand that a thermodynamic quantity like the ground state energy should 
only depend on the numbers A''^, and not on which sector we are working in. 
This is equivalent to requiring that the ground state should not show phase 
separation. Let us therefore calculate the ground state energy in the sector 
where the different species are completely segregated. The A^^^ particles of 
species a occupy a length La, so that J2a = L. For a single species, we 
know that the ground state energy E^a = '^'^^'^9aa^a/i^^a''^a)- Now we 
minimize Eo = J2a Eoa as a function of the La keeping the sum of the La 
fixed; this is equivalent to saying that the pressures exerted by the different 
species must be equal. (In the thermodynamic limit, the contributions to 
the energy coming from the regions separating two species can be neglected). 
This calculation leads to the expression in (pUj). For later use, we note that 
the product rriaNa/La is the same for all a in this sector. In fact, the mass 
density X)a i^aPa must have the same value at all positions and in all sectors 
due to Eqs. ( p!^ - |T3|) , because the pressure must be the same throughout the 
system. However this sum rule does not fix the individual densities pa which 
depend on the sector (see Sec. 4). 

Since particles cannot exchange their positions in the CS model, two 
species are indistinguishable if their masses are equal. (Actually, this is only 
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true for the model on a line, not on a circle. For instance, for two particles on 
a circle, we will show below that the energy spectrum is different for identical 
and non-identical particles). We will assume that all the masses are unequal; 
if some of them are equal, the CS model on a line reduces to one with a 
smaller number of species. 



3. A Multispecies Calogero-Sutherland Model 

A CS model satisfying the relations in Eq. has another attractive fea- 
ture, namely, that its ground state can be solved exactly [jl8|. We first present 



the Hamiltonian for a general CS model defined on a circle of perimeter L. 

Y 2mi f^. ' ' 2mimj sm^[ f (xj - x^) ] 

where the sum over z,j runs over all the N = J2a^a particles. Since the 
Hamiltonian does not distinguish between gij and l — Qij, one needs to specify 
the boundary condition on the allowed wave functions, namely, that ~ \xi — 
Xj\^^i when Xi — Xj —>■ 0. Since the wave functions vanish at coincident points, 
it is clear that there are a number of disjoint sectors which correspond to 
different orderings of distinguishable particles as we go anticlockwise around 
the circle ^ . Under time evolution, a configuration of particles cannot change 
from one sector to another. To illustrate the idea of sectors, consider a system 
with two particles of species 1 and two of species 2. Then there are two 
sectors, namely, (1122) and (1212). 

We will work in a bosonic basis; if all the N particles belong to the same 
species, we demand that 

^(0,Xi,X2, ...,X7v_i) = ^(Xi,X2, ...,XAr_i,L) , (17) 

where < Xi < X2 < ... < a;7v-i < L. (In a fermionic basis, the two sides 
of ( p!7D would differ by the phase (—1)^^^ ). On the other hand, if all the 
particles are distinguishable, we demand that the wave function should be 
periodic in each of the particles separately; namely, 

^(xi, ...,Xj_i,0,a;i+i, ...,X7v) = ^(xi, Xi_i, L, x^+i, x^v) (18) 

whenever < Xj+i <... <X7v<a;i<... < Xj_i < L. We can similarly 
work out the appropriate boundary conditions if some of the particles are 
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distinguishable and some are not. In every case, we see that it is sufficient to 
specify the wave function in the domain < xi < X2 < ... < xat < L; we will 
work in this particular domain from now on. The wave function in all other 
domains can be obtained by using the appropriate boundary conditions. It is 
important to distinguish between the concepts of sectors and domains. For 
instance, in the above example with four particles, the sector (1122) has four 
domains while the sector (1212) has two domains. Under time evolution, 
particles can move from one domain to another. For any and any sector, 
it can be shown that is exactly divisible by the number of domains D. 
Further, a sector can be specified by repeating a pattern of D integers N/D 
times; the D integers are chosen from the set {1, 2, A}. For a single species 
CS model, there is only one sector and one domain. 

For N = 2, there is only one sector; the number of domains is two and 
one if the particles are distinguishable and indistinguishable respectively. 
The energy spectrum can be found exactly by separating the wave function 
into functions of the centre of mass and relative coordinates; however we will 
see that there is an interesting selection rule if the particles are indistinguish- 
able. The wave function is labelled by two integers ni, the centre of mass 
quantum number which runs over all integers, and n2, the relative coordi- 
nate number which must be zero or positive. In terms of the coordinates 
X = {rriiXi + 1712X2) / {mi + 7712) and x = X2 — Xi, we have 

X Tlx 7TX 

= exp ( i27rni— ) | sin (— ) j^^^ Q„,( cos (— ) ) , (19) 

where Qn2 is a polynomial of degree n2', it is even or odd depending on 77,2- 
The energy is 



_ ^ 

-^nin2 — j^2 



2n\ 2 ^1 + 



+ K + 912? -^i^-^ . (20) 



mi + m2 2mim2 



For indistinguishable particles, we set mi = m2 and gi2 = Qu in Eqs. (|T^ 
pO|) ; in addition, we must restrict rii + 77-2 to be even (or odd, depending on 
the choice of basis). To see this, we note that ^'(O, |/) = {—l)^^'^^'^'^{y,L). In 
a bosonic basis, rii + n2 must be even. In a fermionic basis, ni + n2 must be 
odd. Note that the statement that the two particles are bosons (or fermions) 
if 5^11 = (or 1) is independent of the choice of basis. 

For > 3, the ground state of Eq. (|16]) cannot be found exactly in 
general. It is also not clear that the ground state energy will be the same 
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in all sectors. However the situation simplifies dramatically if we impose the 
relations (||). The exact ground state wave function is then given by 



= n I sin [ y {x,-x,) ] 1^- . (21) 

This is the ground state since it has no nodes for Xi ^ Xj. The energy is 
exactly the same as in (p!OD. It is clear that the energy is the same in all 
sectors. 

If we write an eigenstate of Eq. (0) in the form = $\l/o, then $ must 
satisfy the eigenvalue equation 

- [ E ^ + ^ E cot[|(x.-x,)] (m,a. - m,d,) ] $ = (E-E^) $ , 

(22) 

where di = d/dxi. A class of exactly solvable eigenstates is given by the 
excitations of the centre of mass; these are discussed in detail in the next 
paragraph. We have not been able to find any other excited states exactly; 
however we will discuss some approximate low-energy eigenstates in Sec. 6. 

We will now study the centre of mass excitations assuming that the reala- 
tive coordinates problem is in its ground state (^). In our preferred domain 

< xi < X2 < ■■■ < xn L, the centre of mass excitations have the wave 
function 

2T[n „ 

i 

where M = Y^^rrii is the total mass. We now have to find the allowed values of 
n. To do that, we have to consider all the domains and boundary conditions 
together. For a given sector, let us label the different domains by an integer 
d = 1,2,. ..,D, where d = 1 denotes our preferred domain. The functional 
form of the wave function in any domain d is given by 

$d = $1 exp {tOd) , (24) 

where $i is defined in (^) and = 0. The boundary conditions now imply 
that if we go from the domain 1 to a domain 2 by moving only one particle 

1 across the point x = from L~ to O"*", then 

e, - 9^ = 2nn ^ (25) 
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(There is no extra phase arising from the relative coordinates wave function 
(p]) since it remains unchanged if a coordinate Xi changes by L). We can 
then go to a domain 3 by moving another particle across x = 0; if we repeat 
this process D times, we come back to domain 1 with the wave function $i 
times a phase. We can show that the total phase accumulated is 27inD/N. 
Hence nD/N must be an integer. We therefore have the important result 
that the spectrum of centre of mass excitations 

AEn = ( I (26) 

depends on the sector through the number of domains D. We would like 
to emphasize that this discussion of the centre of mass excitations holds for 
a general mulitspecies CS model in which the interaction parameters may 
not satisfy (H); this is because the dynamics of the centre of mass does not 
depend on the gi/s. 

The ground state of the multispecies CS model is exactly solvable if we 
retain (^) but consider the problem on a line in the presence of a simple 
harmonic potential, with the same harmonic frequency u for all the particles 
The Hamiltonian is 

Y f^. 2m.imj [xi - XjY 2 ^ 

The number of different sectors now is A^! /(Ha^a! )• In any sector, the 
ground state wave function is 

= n exp [- 7^ E "^i^i ] ' (28) 

and the energy is 

^(E-«^0' + T^E(l-«-niVa. (29) 

We can find an infinite number of excited states (but not all the excited 
states) as follows. First of all, the Hamiltonian separates into a centre of 
mass problem with coordinate X = rriiXi /M, and a relative coordinate 
problem. The latter further separates into one 'radial' coordinate r where 

= ^E "^i^j {^i-Xjf , (30) 
i<j 
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(this denotes the moment of inertia about the centre of mass) and N — 2 
'angular' coordinates collectively denoted by Qi. The term inside the ex- 
ponential in Eq. (ESI) is J^ifT^i^f = M{X'^ + r^). The Hamiltonian can be 



written as 

2M dX^ 2 

fi^ , N-2 d , 1 2 2 ^ .O.N 

where £ is a differential operator which only acts on functions of Qi. The 
eigenvalues of the centre of mass Hamiltonian are given by huj{ni + 1/2); the 
corresponding wave functions take the form 




*.,(A-) = P„.((/^|;.Y) e.xp ( - ^A- ) , (32) 

where P„j denotes a polynomial of degree rii. The entire difficulty in finding 
the spectrum of ( pl|) lies in finding eigenstates Y{Vti) of C which satisfy the 
boundary conditions at coincident points and the symmetry between particles 
of the same species. For any eigenstate of £ with eigenvalue A(A + — 3), 
where A > 0, we get an infinite tower of eigenvalues of the relative coordinate 
Hamiltonian of the form huj{2n2 + A + (A^ — l)/2). The r-dependence of the 
corresponding wave function takes the form 

^n.(r) = QuA^r^) exp ( - -^r^ ) , (33) 

where is a polynomial of degree it has n2 nodes. The excitations 
corresponding to are called breathing modes. In particular, a set of 
exact eigenstates of (|3l| ) is given by 



^nin. = PnA^^X) Qn,{^r') , (34) 

and the energy is E = Eo + huj{ni + 2n2). 

The Hamiltonians (IT6|) and (ETf) can be written as 



H = Q\Q^ + , (35) 
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where 

1 , inh TT 



Qi = [Pi + — J2 9tj cot - {xi - Xj) ] (36) 



in and 



<3i = 1 — [ Pi + XI ^7 — ~ irriiUXi ] (37) 

ATfli Xi Xj 



in (pTl). The states (|21|) and ( pH]) are annihilated by all the Qi in (|36|) and 
(|37D respectively. 

Incidentally, one can check that some of the assumptions made in Sec. 2 
are indeed valid for the CS model discussed here. Firstly, for two particles 
on a line, with the Hamiltonian being given by ([27| ) but without the simple 
harmonic confinement, we can solve for the scattering problem after going 
to relative coordinates. We find that the scattering phase shift is given by 
'Kgij9{vi — Vj). Secondly, if all the parameters gtj are integers, we can expand 
the ground state wave function (^); it consists of a superposition of plane 
waves with different values of the particle momenta. We then discover that 
all the solutions of Eq. (^ appear in the superposition. There are also terms 
in the superposition which do not correspond to solutions of (^; however 
this is true even for the single species CS model. 



4. Thomas-Fermi Calculation of the One-Particle Densities 

In this section, we will use a generalized Thomas-Fermi (TF) method 
[poll to obtain the one-particle densities Pa{x) in the ground state of the CS 
model in an external harmonic potential. The advantage of using the TF 
method is that one does not need to know the exact wave function; for the 
single species CS model, the semicircle law comes out quite easily. For the 
multispecies model, we will discover that the particle densities depend on the 
sector. 

In the absence of an external potential, Eq. ( |Tl| ) says that 

= irha nib Pb (38) 
ma V 

is independent of both a and x in all sectors. In the presence of a potential, 
the TF method uses Eq. ( pSj ) in a local sense, i.e., ppa and pa become 
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functions of x. The most energetic particle of species a at position x must 
satisfy 



where /i^ is the global Fermi energy for a particles, must be independent 
of X, otherwise one could move an a particle from a position with a higher 
value of /ia to a position with lower /i^. This would give us a different sector 
with a lower energy which is not possible). We thus see that 

must be independent of a. Clearly, the particle densities must vanish outside 
an interval [ 1, where 



X 



2 2/ia 2ha 



Y.^bN,. (41) 



Eq. (1^) only gives us a weighted sum of the particle densities, not the 
individual values of pa- This is because the individual densities on the left 
hand side of P0| ) depend on the ordering of the particles, although the right 
hand side is the same for all orderings. 

If we are given a particular ordering, we can compute the densities as 
follows. Introduce a parameter r, where < r < 1, and functions /a(T) such 
that the first tN particles from the left contain fa{T)Na particles of species a. 
{N is the total number of particles). We assume that /a(T) varies smoothly 
with r. Clearly, /a(0) = 0, /a(l) = 1, and /^(t) > 0. We emphasize that 
the functions fa are fixed once the ordering of the particles (i.e., a sector) is 
specified. Now we want to express r as a function of x. At a; = —Xo and Xo, 
T = and 1 respectively. In the interval x to x + dx or, equivalently, r to 
T + dr, the number of a particles is equal to 

Pa{x) dx = Nafair) dr . (42) 

From Eqs. (^0|-^), we deduce that 

u; r dy {xl - xY^ = TTha rua Nafair) , (43) 

J—Xn 
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so that r is implicitly a function of x. The densities then follow from Eq. 

{^,Pa{x)=NJ',{T)dT/dx. 

It is interesting to observe that Eq. (JI|) is valid semiclassically even in 
the presence of a simple harmonic potential. We simply use Eq. ([T^ ) with 
ha = /2ma + maUp'x^ and let e — after calculating AZ}a(e)- 

We can use the TP method to calculate the one-particle densities even if 
the confining potential is not simple harmonic and the ground state is not 
exactly solvable. \iVa{x) is the one-body potential felt by a particle of species 
a, we can find the condition under which the ground state energy will be the 
same in all sectors. As in Eq. (pQ]), we want 

^ + Va[x) = fia (44) 

with fia independent of x. From (p8D, this implies that (/Xq — K(x))/ma 
should be the same for all species. We conclude that d(ya{x) / rria) / dx (i.e., 
the classical acceleration) should be independent of a. 



5. High Temperature 

The entire analysis so far has been at zero temperature. We now study 
the CS model at finite temperature. Two things are particularly important 
to check. 

(a) Are thermodynamic properties like the equation of state the same in all 
sectors (as is true at zero temperature)? We will argue that the properties 
probably depend on the sector at high temperature. 

(b) Does the model show mutual statistics, if not exactly, then at least ap- 
proximately in some range of temperatures? 

In this section, we consider the equation of state of the system at high 
temperature. The pressure P can then be expanded in terms of the densities 
Pa = Na/ L. We will only be interested in the expansion upto quadratic order, 
i.e., in the second virial coefficients. The grand canonical partition function 
Z has an expansion of the form 

Z = 1 + CaZa + CZaa + E UhZa, + •.• , (45) 

a a. a<b 

where denotes the fugacity of species a, Za is the partition function of 
a single particle of species a on a circle of length L, and Zab = Zba is the 
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partition function of two particles of species a and b. Next, the densities 
satisfy 

(a dlogZ 

We then obtain 

13P = Pa + ^aa pI + PaPb + •.• , (47) 

a a a<b 

where /3 = l/Zc^T, and the second virial coefficients are given by 



hm L ( - 



Aab = hm L ( 1 - ^ ) . (48) 
In the hmit L ^ oo, we can use the Poisson resummation formula 

oo oo 

Y exp ( -7mV^^) = ^ E exp ( - vm^L^) (49) 

n.=— oo n.=— oo 

to show that Za = L/\a plus exponentially smaU terms; = (27r/i^/3/ma)^/2 
is called the thermal wavelength. Using Eq. (^), we can calculate Zaa and 
Zab- We finally obtain 



Aab = {gab - 2) V^a + ■ (50) 

As mentioned earlier, there is only sector for the two-particle problem on a 
circle, but there is an exponentially large number of sectors in the thermody- 
namic limit N ^ 00. The virial coefficients calculated above must therefore 
be some kind of an average over the different sectors, although we do not 
know at this point if all sectors contribute to the average with equal weight. 

Now there is one sector in which we can calculate many properties exactly; 
this is the sector defined at the end of Sec. 2 where the different species are 
segregated, with Na particles of species a occupying a segment of length La. 
We may compute the second virial coefficients in this sector as follows. For a 
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single species, we know from the above equations that the pressure of species 
a is given by 

PPa = Pa + Aaa pl + •.• , (51) 

where pa = Na/La. At equihbrium, the pressures exerted by different species 
must be equal. This enables us to eliminate the La in terms of the total 
length L and the densities pa- We thus get 

f3P = + E Pa + E (^aa + ^bb) Papb + ... . (52) 

a a a<b 

We see that the virial coefficients Aab in this sector do not agree with those in 
Eq. (|50[) since the masses are unequal. Thus the equation of state seems to 
depend on the sector. Since the second virial coefficients should be related to 
the statistics parameters |P, |l^, we conclude that the model cannot realize 



mutual statistics in all sectors at high temperature. 

We now observe that the particle density J2a Pa is independent of both 
the position x and the sector at T ^ cxd, while the mass density J2a^aPa is 
independent of the position and sector at T — > (cf. Sec. 2). This implies 
that the individual densities pa must vary with the temperature in at least 
some sectors. Such a variation may seem surprising in a model which is scale 
invariant, in the sense that the kinetic terms and the two-body interactions 
transform in the same way under an uniform scaling of all the coordinates. 
However the temperature-dependence might be understandable if the model 
is found to be non-integrable. We have therefore checked for integrability, 
both classically and quantum mechanically, using the methods given in Refs. 
[|, H; the model appears to be non-integrable. Hence the asymptotic Bethe 
ansatz is not expected to be valid for the excited states, except for the centre 
of mass excitations in which all the velocities are shifted by the same amount 
from their ground state values. 

We recall that in the single species CS model, the exclusion statistics is 
ideal. This happens because the scattering phase shift is independent of the 
momenta of the two particles and the model is integrable 0, ^. The second 
property implies that an A^-particle scattering phase shift is given by a sum 
of two-particle phase shifts. 

Although the CS model does not seem to exhibit mutual statistics in all 
sectors at high temperature, it may do so in an approximate way at low 
temperature. In the next section, we study the low-energy excitations in 
some detail. 
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6. Low Temperature 



A general formalism for obtaining the distribution functions for a system 



with mutual statistics is given in Ref. |]14|. It is therefore interesting to 



investigate, (a) whether the low-energy excitations are described by one or 



more c = 1 Gaussian conformal field theories |jT^, and (b) whether the 
equation of state of the CS model at low temperature agrees with that of a 
system with mutual statistics. We will provide a complete answer to (a), but 
we are unable to answer question (b) by analytical methods. 

The elementary excitations in a Gaussian theory satisfy the dispersion 
relation 

= vf\^p\ (53) 

where AE and Ap respectively denote the energy and momentum with re- 
spect to the ground state, and vp is the Fermi velocity. From Eq. ([T^), 
vp = nhaM/L, where M = '^^i^i is the total mass; vp remains finite in 
the thermodynamic limit M, L — >■ oo. (Actually, Eq. ( |5^ ) describes the ex- 
citations of a Gaussian theory based on an uncompactified scalar field. For a 
compactified scalar field, there are some additional terms in the energy which 
will be discussed in Eq. (|66D below). 

We can now find a right-moving elementary excitation approximately by 
using the Sutherland ansatz for the wave function 0. Let us first define 

Zi = exp{i—Xi) . (54) 
ij 

Following Eq. (^2|), we write \E'i = $i\l'o where 

i 

This has Ap = 27Th/L since the momentum operator is p = J^iPi- We can 
show that the expectation value 

{^i\H-Eo\^i) Ap 
= ^-^-^ = vpAp + 0(— ) . (56) 

The corrections of order Ap/L arise due to the approximate nature of the 
wave function, and they go to zero as L,N —>■ oo; we then recover AE/Ap = 
Vp. Similarly, a left-moving elementary excitation with Ap = —2nh/L is 
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given by = <l'_i\&o where $_i = It is important to note that the 
leading term in AE / Ap is independent of the sector. 

The first term in AE/Ap in Eq. (|5B[), which survives as L ^ oo, arises 
entirely from the /mear differential operators in Eq. (^2]). On the other hand, 
the corrections of order 1/L arise from both the linear and the quadratic 
differential operators in (^21). It is very useful to know that the leading term 
in AE arises from a linear operator; it implies that the energy eigenvalues 
will add up when we multiply two eigenstates together. Thus a family of 
composite excitations with Ap = 2TTh{ni — n2)/L is given by 

^nun, = $r -^-'l , (57) 

where ni,n2 > 0; this state has AE = VF^iih^ni + n2)/L. We require that 
ni,n2 << N so that this is a low-energy excitation and the corrections of 
order Ui/L vanish as L — oo. 

Following Ref . , we can show that the general excitation is specified by 
a set of non-zero integers (ni, n2, n^); they may be positive or negative. 
(For this to describe a low-energy excitation, we again require that each of 
the magnitudes \ni\ as well as k should be much less than A^). The energy 
and momentum of such a state are given by 



AE = Vf ^ y 

Ap = -—J^m. (58) 

^ i 

If all the integers rij are positive, the ansatz wave function ^ni,n2,...,nk contains 
terms like {miZi)"'^ (771222)"''^ ■■■{mkZk)"''' plus permutations, plus other terms 
which are obtained by squeezing (see below). If an integer Ui is negative, 
{rriiZi)^' is replaced by {miZ*)~'^\ The wave function $ni,n2,.--,nfc is obtained 
by multiplying together the ansatz wave functions of a number of 'elementary' 
states; an elementary state is one in which all the integers rti are equal to 1, 
or all of them are equal to —1. 

Let us first examine the wave function for an elementary state, say, (1, 1). 
This is given by \E'i 1 = $11 where 

$1,1 = y f{xi - Xj) rriimjZiZj . (59) 
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The factor f{xi — xj) in ( |59D requires some explanation. If the particles 
i and j belong to the same species, we set / = 1. But if they belong to 
different species, we require that / should go to zero sufficiently fast if i and 
j approach each other, and should go to 1 if their separation is much larger 
than the inverse density L/N. For instance, we could define the function as 

The reason for introducing such a function in (|59|) is that if we always set 
/ = 1, the linear operator in Eq. (|22| ) would act on ( p9D to produce a 
function which diverges when Xi — >■ Xj; the expectation value {"^IHI"^) would 
then be ill-defined if Qij is not large enough. (We would like to emphasize 
that the function / is only required because we are working with approximate 
wave functions; such a short-distance cutoff is not present in the exact wave 
functions). Since / only differs from 1 in certain microscopic regions of 
space, its presence can only change the expectation value of AE/Ap by an 
amount which vanishes as L —>■ oo. We can now verify that the energy and 
momentum of the state (^) is given by (^8]). 

We can write down the wave functions of other elementary states like 
(1, 1, 1) in a similar way. The wave functions of left-moving states such as 
(— 1, —1) are obtained by complex cojugating Finally, composite wave 
functions such as $2,1,-1 can be written as a product of $1^1, $1 and The 
general rule for writing $ni,n2,...,nfc as a product of elementary wave functions 
can be expressed inductively as follows. Suppose that the numerically largest 
positive (or negative) integer appears / times in the set (ni, n2, n^). For 
instance, let ni = n2 = ... = n/ be the largest positive integers. Then 

^ni,n2,...,ni,ni+i,...,nk = '^'ni-l,n2-l,...,n;-l,n;+i,...,nfc '^'l,l,...,l ; (61) 

where the last $ has / I's. This factorization can be repeated till we only get 
elementary wave functions on the right hand side of Eq. (pT|). Ifrii = ... = ni 
is negative, we replace n — 1 by rii + l, for 1 < z < Z, and by $_i 

on the right hand side of (0). An important point to note in all this is that 
a factor of comes along with each Zi or z*. 

Finally, we observe that $2,1,-1 contains within it other wave functions 
such as $1,1,1,-1, $1,1 and $2, because the sets (1, 1, 1, —1), (1, 1) and (2) can 
be obtained form the set (2, 1, —1) by a process called squeezing 0. We say 
that a set {n-} is obtained from another set {rij} by squeezing if, either 



20 



(a) an integer rii in the second set is equal to the sum of two integers n'^ and 
n'f^ in the first set, where rzj, n'^ and n'^ all have the same sign, or 

(b) there are two integers nj,n; in the second set, and upto two integers n'^ 
and n'^j, in the first set, which satisfy ni + rii = n'j + n'j. and rii > n'j, n'f^ > rii. 
(One or both of the integers n'- and n'^. may be zero, in which case it is 
omitted from the set {n[}). 

Following (|58|) , we can think of the low-energy excitations in terms of a 
system of 'bosons' whose allowed energy levels are given by 

en = vp -J- \n\ , (62) 

where n can be any non-zero integer. These bosons have zero chemical po- 
tential. The partition function is therefore 



exp( 



r 1-2 

-f3Eo) [ n ( 1 - exp(-/5e„) ) ] . (63) 



We can now calculate the leading terms in the specific heat and pressure 
at low temperature. The mean energy is given hy E = — {din Z/dlS)^, the 
specific heat is Cy = {dE/dT)^, and the pressure is P = ksT^dln Z/dL)T- 
The leading terms in the specific heat and pressure are therefore given by 



Cv 7Tk%T 



V 



F '"^B 



37rha 3hvF 



(64) 



These leading terms are independent of the sector. However we are unable to 
obtain the next term of order in Cy or order in P, since that requires 
a knowledge of the 1/L terms in the energy; we have ignored those terms 
because they do not seem to be analytically computable. The importance of 
the next term in Cy or P is that it contains information about the statistics 
parameters |T5 . 



Thus the low-energy excitations match those of a c = 1 Gaussian con- 
formal field theory, rather than c = A (as would be the case if all the off- 
diagonal elements gab were zero). This raises the following question. Are the 
low temperature thermodynamic properties of the model the same as those 
of a single species CS model with some 'effective' particle density p, mass m. 
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and interaction parameter g7 On comparing the specific heat and pressure 



in (|64D with those of a single species model, we see that the three effective 



parameters must satisfy the two relations 

nhgp 



m 
_9_ 



a . (65) 



Clearly, we need a third low-energy property (for instance, the next term 
in Cv or P) to fix the individual values of p, m and, most importantly, 
the exclusion statistics parameter g. As mentioned earlier, we also have to 
determine whether or not these values depend on the sector. 

Before concluding this section, we note that the centre of mass excitations 
discussed around Eq. (^) can have energies which are much lower than 
the ones discussed above if the number of domains D in a given sector is 
comparable to N . If the number of domains D « N (as in a single species 
CS model where D = 1), the centre of mass energies are of the same order 
as the ones discussed above. In any case, the centre of mass constitutes only 
a single degree of freedom, so it has no effect on thermodynamic properties. 

For the sake of completeness, we should mention that the low-energy ex- 
citations of a single species CS model have some terms in the expressions for 
energy and pressure in addition to the ones given in (pH]). These are char- 
acterized by two integers, AK equal to the number of particles transferred 
from the left Fermi momentum —pp to the right Fermi momentum pp, and 
AA^ equal to the change in the number of particles. The contributions of 
these to the energy and momentum take the form []T1| 



AE = i^AN + vp -— ^ + ^{AN) 



L 



\2 



9 4' 



Ap = 2-nhp AK + ARAN . (66) 

We observe that AK simply corresponds to a centre of mass excitation in 
which all the particle momenta are shifted by the same amount 2tt%AK / L; in 
the notation of Eq. (^6l), AK = n/N. Eq. (|66|) is essential in order to identify 
the radius of compactification of the scalar field in the Gaussian theory as 
R = In the multispecies model, we can certainly calculate terms 

similar to (^) using Eqs. (^) and (|10[); however we will not exhibit them 
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here since these terms have no bearing on the thermodynamics poperties 
which were calculated above with fixed particle numbers Na. 



7. A Multispecies Anyon Model in Two Dimensions 



In this section, we would like to point out that relations similar to (|]) lead 
to some special properties for mutual statistics in two dimensions. Consider 
a collection of A species of anyons with mutual statistics parameters gab = Qba 
and charges [||. Let us first assume that all the have the same sign 
and all the Qab are positive. We now apply an uniform magnetic field in a 
direction normal to the plane. At zero temperature, all the particles lie in 
the lowest Landau level (LLL). Let Bqa > 0. Then in the symmetric gauge 
{Ax, Ay) = {—By/2, Bx/2), all the wave functions in the LLL have the form 



f{zi) exp 



B 



J2 li'^i' 



(67) 



where Zi = Xi + ii/i and rf = ZiZ*; c is the speed of light. (Note that the 
particle masses play no role in the LLL). When anyon i is taken around a 
closed loop surrounding anyon j, the wave function should pick up a phase 
exp{i27rgij). The simplest wave function of this kind is given by 



f{Zi 



n (^^ 

Kj 



(68) 



This is multivalued if any of the gij is not equal to or 1. Other wave 
functions in the LLL are obtained by multiplying (^) by any polynomial in 
the Zi which is symmetric between anyons of the same species. We will only 
consider the wave function (|67| - |68|) henceforth. Amongst all the LLL wave 
functions, this has the lowest angular momentum, i.e.. 



L = { Zi 



d_ 

dzi 



- z, 



d_ 

dz* 



) 



(69) 



is equal to fiY,i<j Qij- The particles with this wave function lie closest to the 
origin Zi = 0. 

The probability density is given by 



1^1 



n 

t<j 



exp 



B 
2ch 



(70) 
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We now note that the expression (|70D may sometimes be vahd even in a 
situation in which the wave function is not given by (|67D; in that case, qi 
in (^) will denote the magnitude of the charge of species i, rather than the 
charge itself. This is in fact true for the example given at the end of this 
section. 

We will calculate the one-particle densities Pa{r) using a plasma analogy 
[p]] . Since ([70|) is circularly symmetric, pa can only be a function of the 
magnitude r. Now — ln|\E'p can be interpreted as the potential energy of 
a system of point charges interacting with each other through a repulsive 
Coulomb potential (which is logarithmic in two dimensions) and with an 
uniform attractive background charge density equal to —B/nch. A particle 
of species a which is at a distance r from the origin will be at equilibrium if 

E 2g., ^ = -f-lar, (71) 

where ^ 

Q,{r) = 2n [ r'dr' p,{r') (72) 





is the total charge of species b inside the circle of radius r. Eqs. (|7l]-0) 
imply that 

vr E 9ab Pb{r) = ^ 9a ■ (73) 

This is a system of A equations if all the densities pa are non-zero at r; if 
some of them are zero, we have to drop the corresponding equations from 

We now demand that the system should satisfy the following require- 
ments. Firstly, since the right hand side of (|7BD does not depend on r, we 
would like the particle distribution to be homogeneous, that is, pa should 
be independent of r. Secondly, if the numbers of particles Na are large but 



finite, we demand that there should be a single radius R» ^ch/Bqa such 
that all the densities pa are non-zero inside R and zero outside R. These two 
requirements rule out the possibility of phase separation, namely, a situation 
in which some, but not all, of the densities vanish in some region. We then 
have 

BR^ 

E 9ab Nk = Qa . (74) 
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Finally, we demand that the numbers Na be independent variables; for in- 
stance, A^i should not get fixed if N2,N3, ...,Na are specified. We can show 
that these three requirements will be satisfied if and only if all the equations 



in (^) are linearly dependent on just one of them, i.e., if 

gab = a QaQb (75) 

for both a = b and a ^ b. 

It is interesting to note that the relations in ( [75|) arise quite naturally 
in the fractional quantum Hall effect. For instance, consider the situation 
when the filling fraction is close to but slightly less than 1/m, where m is 
an odd integer like 3 or 5. In the first hierarchy, the ground state contains 
two species of particles with opposite charges, namely, the electron (species 
1) and the quasihole (species 2). Laughlin's wave function then tells us that 
gii = m, gi2 = I, (?22 = and the ratio of the magnitudes of the charges 



is q'2/Q'i = 1/"^ hence (^) is satisfied. Ref. [§] points out the interesting 
fact a fractional quantum Hall system can be described by a single Chern- 
Simons gauge field only if Eqs. (^) are satisfied. 



8. Discussion 



To conclude, we have argued that certain relations between the interac- 
tion parameters and masses are necessary for the ground state of a one- 
dimensional scale invariant multispecies model to satisfy the asymptotic 
Bethe ansatz. We have shown that a CS model with those relations has 
an exactly solvable ground state, and that the low-energy excitations are the 
same as those of a Gaussian conformal field theory. However we have not 
proved that the same relations are sufficient to give rise to either mutual 
statistics, or to a single species exclusion statistics. To address that ques- 
tion, we need to understand the spectrum of low-energy excitations more 
accurately. 
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